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Abstract 

The electromagnetic contribution to the isomultiplet mass split- 
tings of heavy mesons is reanalyzed within the framework of the heavy 
mass expansion. It is shown that the leading term in the expansion is 
given to a good approximation by the elastic term. 1/mg-corrections 
can only be estimated, the main source of uncertainty now being in- 
elastic contributions. The l/m<g-corrections to the elastic term turn 
out to be relatively small in both D and B pseudoscalar mesons. 
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1 Introduction 

The measurement of mass splittings in heavy pseudoscalar and vector mesons 
has attained a substantial degree of precision in the D-system |], ||] and in 
the B-system is expected to improve over the presently available data || . In 
particular, the isospin splitting in D and D* mesons is now known within an 
error of ±0.2 MeV f2|, in B mesons the error is larger and about ±1 MeV 
[|3| and, in B* mesons its experimental determination is an important task 
still to be accomplished. The theoretical understanding of isospin splittings 
is less satisfactory. Neither the electromagnetic nor the quark mass contri- 
butions can be determined to a similar degree of precision. In both cases, 
the task boils down to a strong interaction problem where the extraction of 
a result cannot be at present achieved without assumptions and, to some 
extent, modeling f|, ||, ||, |7|, j|, ||]. Surprisingly good agreement with data 
was found || by an analysis in the limit of infinite heavy quark mass using 
a linear interpolation to determine the contribution A m M to the splitting 
caused by the mass difference of the u and d quarks, and the elastic approxi- 
mation to the electromagnetic contribution A 7 M with a VMD model for the 
elastic form factors associated with the light quark components of the elec- 
tromagnetic current. There are, however, unanswered questions concerning 
each of the assumptions made. In the present work, we reconsider A 7 M in 
the light of the heavy mass expansion. We consider models beyond VMD and 
give an estimate of corrections of order I/ttlq. Some conclusions of general 
character are drawn from the analysis. 

2 Electromagnetic mass splittings: pseu- 
doscalar mesons 



The QED contribution to the mass of a hadron H is given to order a by: 
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where P 2 = M^, — iD M „(g) is the photon propagator in an arbitrary gauge, 
the tensor T^ the covariant forward amplitude for virtual Compton scatter- 
ing off H, given by 



T%(q, P)=i d'x e'** (H, P,<r\T f(x)f(0) | H, P, a) 



(2) 
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where a denotes the polarization of the hadron, and S^M^\S) is the coun- 
terterm needed to render 5 7 M# U.V. finite. The customary covariant normal- 
ization of states is used. In this work, we focus on heavy mesons containing 
a light quark (u, d) and a heavy antiquark (c, b). H and H* will denote re- 
spectively the pseudoscalar and vector mesons. The electromagnetic current 
pertaining to our analysis is conveniently decomposed as follows: 

f = JS + \-% + eQ-% 

4 = \{uru-d^d) (3) 
Jy = \ (u^u + drfd) 

J" Q = QYQ 

where Q denotes the heavy quark with charge cq. 

Since QCD plus QED is a renormalizable and parity conserving theory, 
the counterterm £ 7 M|^'(A) can only involve the matrix elements of even 
parity gauge invariant operators of dimension four or less [TO]. At order a, 
the only such operators are G^G^ {G^ u is the QCD field strength), and 
QQ. The quark bilinear operators must appear multiplied by a coefficient 
proportional to the respective quark mass as demanded by a consistent chiral 
limit. In this limit, the difference of mass shifts relevant to this work 

A 7 Mjj = <P Mu d - 5m Hu (4) 

becomes U.V. finite. This results from the cancellation of the gluon- and 
heavy quark- operator contributions in the difference. Throughout this 
work, we take the chiral limit; this introduces an error in A 7 M# of order 
am q / A.QCD- The sj^s component of the electromagnetic current drops off 
when considering (4). For this reason it was not included in (3). 

In this section we discuss the case of pseudoscalar mesons. Following 
Cottingham fLT |, one can perform in (1) an Euclidean rotation in the v = ^ L 



complex plane and express: 



where Q 2 = —q 2 , and: 
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It is convenient to separate T M1/ into two pieces, one referring only to 
the light quark clouds and the other to the light quark cloud and the heavy 
quark: 



rr + r 



13 



(7) 



d'x e^ x ((H d | T(Jf (0) J»(x) + 4(0) Jf (x)) | H d ) - d 



u 



where % — A, B, and J\ = \ J Y and J V B = cqJq. 

The lack of experimental access to the absorptive part of T^, which 
would permit the reconstruction of T^ v , compels us to resort to some ap- 
proximations. The main approximation will consist in taking only the elastic 
terms in the absorptive part of T^ u . In the framework of the heavy quark 
expansion, we will show that inelastic terms are of leading order [0(1)] for 
i = A and of order 1/mg for i=B. 

The elastic contributions are conveniently studied in the large mq limit, 
where the pseudoscalar and vector mesons belong to a multiplet under the 



heavy quark symmetry [12]. It is therefore appropriate to keep as elastic 



terms in the absorptive part of not only the strictly elastic terms, but also 
those inelastic terms which involve a Ml transition to H* in the intermediate 
state. We can then write the elastic absorptive part of T-^: 
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with the matrix elements 



(Hp,p'\jm\H a ,p) 

(H f) ,P'\Jg(0)\H a ,P) 

(H;,e,P'\Jt(0)\H a ,P) 

(H*p,e,P'\ 4(0) | H a ,P) 

1 
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where e u is the polarization four- vector of H*. In the large ijiq limit one has 
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where £ is the universal I-W form factor, and v denotes the four-velocity of 
the meson. The sign convention for gg given by this equation implies that 
sgnF y =sgns(y. 

In the customary tensorial decomposition of the hadronic tensor, (8) and 
(9) give the following expressions for the elastic structure functions: 

Wtr^P) = -^M 2 H 5(2M H u-Q 2 )e((P + q ) ) gi (Q 2 )g 3 (Q 2 )(u 2 + Q 2 ) 
W$r(Q,P) = -^M 2 H 5{2M H u-Q 2 )6{{P + q f) (11) 



x 
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At fixed Q 2 > the dispersion integral over the variable v gives: 
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Using that: 
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Tftq, v) = -3 T X (Q 2 , v) + T 2 (Q 2 , v){\ + — ) (13) 
replacing in (5), and performing the integration over dv we finally obtain: 
A^ cl M H = Al cl M H + Al cl M H 
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To leading order in the heavy quark expansion the result becomes physically 
very simple: 



A^ l M H = -— dQ (-F Y (Q 2 )F 3 (Q 2 ) 

IT Jo \2 

- e Q F 3 (Q 2 )-^-g Y (Q 2 )g 3 (Q 2 )\ 



(15) 



where the form factors are taken in the mg — > oo limit. The first term is 
the difference between the Coulomb self-energies of the light quark clouds 
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with the u and d flavor quantum numbers, the second term is the difference 
between the Coulomb interaction of these clouds with the heavy quark, and 
the third term is the result of Ml transitions H <-> H* / y mediated by the 
light quark component of the e.m. current. Since the relevant contributions 
to the integral in (15) stem from values of Q 2 of the order of light hadronic 
masses, one can neglect the Q 2 -dependence in Fq(Q 2 ) which amounts to an 
error of order wig 2 . Since Qq is already of order 1/mg, it does not appear in 
(15). 

There are two types of 1/mg-corrections, kinematic ones which are ob- 
tained by expanding the terms in the square brackets in (14), and dynamical 
ones which reside in the heavy mass dependence of the form factors associ- 
ated with the light quark currents. The latter are hard to determine; charge 
conservation assures us that the corrections vanish at Q 2 = and can only 
affect the charge radius. Kinematic corrections due to the mass difference 
Mh* — Mh oc 1/toq are of higher order. One thus obtains the following 
1 / mQ-corrections: 



A\ el Mn 1 = "-^r I dQ [Q F Y (Q 2 ) F 3 (Q 2 ) - ^-Q A g Y (Q 2 ) g 3 (Q 

+ lQ 2 9y(Q 2 )^ r g 3 (Q 2 )} (16) 

Mh 
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Al eL M H = -3e Q — — / dQ (Q F 3 (Q*) - =-Q 2 g 3 (<p 
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We have assumed here that the mQ-dependence of F Y (Q 2 ) and F 3 (Q 2 ) are 
the same. A similar assumption applies to the Ml form factors. In the limit 
where the light quark is taken as a constituent quark, the term proportional 
to gz{Q 2 ) in A B el Mn is identified with the difference between hyperfine split- 
ting. This clarifies our inclusion of the Ml transition terms and also the sign 
of the contribution. 

It is important to notice that the expansion in powers of I/tuq we have 
implemented corresponds to that of the HQET (heavy quark effective the- 
ory). The integration over Q must then be cut off at a scale A << mg, and 
the contributions which arise from the domain Q > A must be included as 
local counter-terms of order 1/mg. These terms are, moreover, responsible 
for lifting the cut off dependence of the results. A list of the possible counter- 
terms one can add has been given in ref 0. The size of the coefficients in front 
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of these terms is not known a priori. In the case of the elastic contributions 
they could be determined once all the form factors are known above the cut 
off. Our ignorance about the counter-terms is manifested, therefore, in the 
cut off dependence of our results, which depends on the asymptotic behavior 
of the intergrands in (16), and which turns out to be more pronounced for 
i = B. 

There are other terms of order 1/toq which have been disregarded; they 
are the tadpole terms resulting from seagull type contributions to T^, and 
which must be present to maintain gauge invariance. In particular, their 
strength is fixed at Q 2 = in terms of the electric charge of the meson. 
This fact has the particular consequence that at order I/tjiq the tadpole 
contributions to A 7 M# and A 7 M#» are equal. 

In conclusion, we are able to determine the 1/mQ-corrections of kinemat- 
ical origin, up to the mentioned cut off dependence, and those originating in 
hyperfine type interactions. Those residing in the form factors will in any 
event be buried in our lack of knowledge of the form factors themselves, and 
those of tadpole origin will be disregarded. 

3 Results 

For the sake of providing a rough quantitative estimate of the mass splitting 
we now resort to a model for the form factors. In the case of the proton 
and the pion, the high Q 2 behaviour of the form factors seems to pervade 
the behavior at lower values of Q 2 , giving rise to the respective dipole and 
monopole shapes. For heavy mesons the situation is more involved. Within 
the heavy quark expansion, the asymptotic behavior corresponds to Aq CD <C 
Q 2 <C m 2 Q. To leading order in a S) the asymptotic behavior of the charge 
and Ml form factors coincide with those of the relativistic hydrogen atom 
with a = |a s : 

F(Q 2 ) oc Q g(Q 2 ) oc Q~ 5 ^ 

a s = a s (Q 2 ) (17) 

The exponent has the following values at some illustrative points: 5(0) = 3, 
5(0.3) = 2.8, 5(0.5) = 2.5, 5(0.6) = 2.2, (5(0.64) = 2. For those values of a s 
where one can trust this estimate the form factors at large Q 2 fall off faster 
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than for the pion. Since the asymptotic behavior is not definitely known, we 
choose the form factors to have the following form: 

3F Y (Q 2 ) = 2F 3 (Q 2 ) = 




3g Y (Q 2 )=2g 3 (Q 2 ) = 

The choice of m p as the relevant scale is natural, and likely to be a reasonably 
good approximation. The constant (3 determines the rate of the transition 
H* — > Hj. We will choose here a value suggested by quark models [f|: (3 ~ 
3 GeV^ 1 . Here k will be varied between 1, which corresponds to the vector 
meson dominance model ||, and 1.5, which asymptotically corresponds to 
a loosely-bound light quark. Neglecting the 1/rriQ corrections to the form 
factors one obtains: 




Al el M n ■ u,— -{-^^-^ (19) 

r(«-i) 



+ ttt" h + P™ 



'' 3 T(k + 1 




+ . . . x 

For the chosen range of values for k, only A^ el M# shows non-analytic con- 
tributions when setting A ~ M H . 



Numerical results for three different choices of k are shown in the table. 
In all cases A\ el Mx is only a fraction of Ag el M#. For (3 = (i.e., in the 
strictly elastic approximation), the considered 1/mg-corrections to A^ el M# 
are small: 0.1 MeV for D mesons and 0.03 MeV for B mesons, while the 
corrections to A^Mh are larger: -0.7 to -0.5 MeV for D mesons and -0.3 
to -0.1 MeV for B mesons, and tend to suppress the splitting. For f3 = 3 
GeV -1 , the additional contributions to A^ el M# and A]j el M# due to the 
Ml transitions are of order mPq and 1/rriQ respectively. They lead to a 
suppression of A^ el M# by a factor 1/3 to 1/2; this, however, affects the 
overall mass splitting by less than 0.35 MeV. The corrections to A]j el M#, 
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identified before to be of hyperfine nature, have the sign opposite to that of 
the kinematical 1/mQ-corrections, and they are about 0.8 to 0.5 MeV for D 
mesons and about -0.3 to -0.2 MeV for B mesons. The cut off dependence is a 
function of (3: A 7 el M# increases (decreases) with increasing A if f3 = {(3 = 
3 GeV -1 ), and according with (19) the sensitivity diminishes as k increases. 
We estimate that the cut off dependence affects the elastic term of D and B 
mesons by ±0.2 Mev and ±0.1 Mev respectively. The dependence on (3 is 
important in D mesons, due to the fact that the dependences of A\ el Mu and 
A^Md add up, while the contrary occurs in B mesons. A 7 el Mo increases 
by about 1 MeV as (3 increases from to 3 GeV -1 , and A 7el Me varies 
by at most ±0.1 MeV in this interval. Overall, for the chosen value of /3, 
1/mQ-corrections are small even for D mesons, and most of the uncertainty 
in A^'Md^ resides in our lack of theoretical knowledge of the form factors. 

At this point, it is appropriate to briefly address the problem of deter- 
mining the total isospin splittings. A m M# is determined from the mass 
difference Mh s — Mn d , which after the inclusion of chiral corrections reads 
& 

M Hs - M Hd = C H (m s - m d ) - (2 M\ - Mj) + 0(m 2 s ) (20) 

where g determines the H* — > Hn amplitude, F =93 MeV is the pion decay 
constant in the chiral limit, and Ch is an unknown constant. Using the ratio 
[7j R = {m s — rh)/ (m^ — m u ) ~ 43, with m the average of the u and d quark 
masses, one obtains: 

A m M H = {M Ha - M Hd + (2 Ml - Mj) ± 0(m 2 s )) (21) 

Thus, the leading chiral corrections, which are non- analytic in m s , tend to 
increase the magnitude of A m M#. Clearly, (21) makes sense only as far as g 
is small enough; the corrections are large for g > 0.3 (for determinations of g 
in the quark model and in QCD sum rules see refs. |16| and for constraints 
see refs. |]17|]). Otherwise, the unknown 0(m? s ) terms must also be included. 
We conclude that there is a large degree of theoretical uncertainty in A m M#. 

For g = and the ratio R ~ 43 one obtains A m M D ~ 2.3 MeV. The 
positivity of the non-analytic term and the experimental result AMh = 4.8 
MeV M imply that A 7 M D < 2.5 MeV. 
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4 Inelastic contributions 



The main omission of the analysis presented here is that of inelastic contri- 
butions. We only kept those involving H* in the intermediate state. The 
inelastic contributions to A\Mh appear at leading order in the heavy quark 
expansion. Since A^ el M# is substantially smaller than A^ c1 Mh, one might 
hope that they will lead to numerically small contributions to the mass shift. 
As already mentioned, those considered amount to less than 0.35 MeV. As a 
curiosity, in the case of the p-n electromagnetic mass difference it has been 



noticed long ago [14| that inelastic terms give only a modest contribution. 
On the other hand, for A^Mh we can prove the following statement: the 
inelastic contributions to A^Mh are of order 1/uiq. The proof is as fol- 
lows: inelastic terms in the absorptive part of involve products of matrix 
elements of the form: 

(H,P | j£(0) | P';X){P';X | J$(0) | H,P) (22) 

where the intermediate state involves a heavy hadron with momentum P' 
and light hadrons denoted by A carrying momentum Px- Let us analyze the 
matrix elements separately. For (P';X \ Jq(0) \ H,P) we use Bjorken's sum 



rule 15 



1 = ^(1 + v.v') \£(v.v')\ 2 -l / dc^c.r.r') (23) 







where u ine \ is given by: 



wUe, v.v') = e7 (o d ^ P> &*f?{P -q-P'-Px) \ (n, P'; X \ J% \ H, 

n;X J V Z7F / Z ^0 

+ Z — graph 

where n denotes a heavy hadron and the prime on the sum implies that for 
A =| 0) one must take n ^ H. Since each term on the RHS is positive, and 
the form factor satisfies the normalization condition |£(1)| = 1, taken near 
the point of zero recoil, Eq. (23) implies the Cabibbo-Radicati type sum rule 



de u inel {e, v.v') = f - + 2 |f |' (1) j (1 - v.v') (25) 

which, supplemented with the assumption that £' is non-singular at zero 
recoil and mild requirements of continuity for coWi, implies that 

[n, P'; X | J£(0) \H,P)\<x (v.v' -If r] > \ (26) 
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The behavior of {H, P | J3 (0) | P'\ X) for v ^ v' is determined as follows: 
in order to change the velocity of the heavy quark by a finite amount the 
large space-like momentum transfer q must flow as shown in the figure. The 
magnitude of the momentum is Q 2 ~ Am 2 Q(v.v' — 1) in the cases where 
P\ -C ttiq, which are indeed the case in our problem. As illustrated in the 
figure, the flow of the momentum q gives rise to a suppression factor of the 
order of 1/Q 3 stemming from the light quark and gluon propagators. This 
clearly holds for all possible diagrams contributing to the matrix element. 
Thus: 

I (H, P I j£(0) I P'; X) I cc 3 ) (27) 

nTq [v.v — I) 6 '* 

which is valid as far as uiq (v .v'— l) 1 / 2 >> Aq C d- Thus, according to (26) and 
(27), in this regime the product (22) becomes proportional to ttiq 3 (v.v'—1)~ x . 
When rriQ {v.v' — I) 1 / 2 ~ ^qcd the LHS of (27) is not suppressed, and (22) 
becomes proportional to ttiq 1 since (v.v' — 1) 1//2 oc rrig 1 in (26). This latter 
case corresponds precisely to the inelastic term, which we included as elastic, 
with H* in the intermediate state (hyperfine term). We therefore conclude 
that Tgf (q, v) is suppressed with respect to the strictly elastic one by a 
factor tuq 3 (v.v' — which upon replacement in (5) leads to the stated 
1/m.Q suppression for A]^ meI M#. 

5 Vector mesons 

Finally, a few comments on vector mesons. Their isospin mass splittings are 
identical to the ones of pseudoscalar mesons at leading order in 1/mg. They 
differ at order 1/mg, where a number of terms can be identified as possi- 
ble contributions || . In particular, for the electromagnetic component there 
are the following terms: the hyperfine term, which breaks the heavy quark 
symmetry and contributes a factor (-1/3) times the hyperfine term for the 
pseudoscalars (the relative factor is the same as in quark models); charge 
form factors F^y f° r pseudoscalar and vector mesons differ at this order (we 
expect this effect to be smaller than the hyperfine term); 1/mg corrections 
to the elastic dipole form factor of the light quark components of the elec- 
tromagnetic current between vector meson states which lead to corrections 
to AY l M H * (expected to be relatively small). The electric-quadrupole tran- 
sitions give corrections of order l/mn, and, as mentioned before, the tadpole 
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contributions are the same as for the pseudoscalars up to corrections of order 

1 / m Q- 

In D mesons the data is summarized by |2|]: 

(a) M D . - M Da = 141.5 ± 1.9 MeV 
M D «o-M D o = 142.12 ±0.07 MeV 

M D *+-M D + = 140.64 ±0.1 MeV (28) 

(b) M D *+-M D , = 3.32 ±0.1 MeV 

M D + - M D o = 4.8 ± O.lMeV 

The results (a) show that the QCD hyperfme splitting is remarkably SU(3) 
symmetric, as emphasized in ||: there is a change of less than 2% as the 
quark masses vary from m u ^ to m s . Unless this is an accident resulting 
from the fact that non-linear terms in m s are large, one concludes that M 
A m M D = A m M D , to a good degree of precision. Thus, AM D , - AM D ~ 
A 7 Mo. — A? Mo ~ —1.48 ± 0.15 MeV, which suggests a hyperfme term in 
A 7el M£> of about 1 MeV, a value which is not very different from the ones 
we obtain, which for (3 = 3 GeV -1 vary between 0.5 to 0.8 MeV. Clearly, the 
possibility of isolating the hyperfme contribution provides a useful constraint 
for further theoretical understanding of I/tuq corrections. 

6 Conclusions 

Electromagnetic contributions to mass splittings in D and B mesons are con- 
veniently analyzed in the framework of the large mass expansion. Leading 
order terms are determined, up to inelastic corrections to A\Mh, by elastic 
contributions which are given in terms of the elastic form factors associated 
with the light quark components of the electromagnetic current. A precise 
theoretical determination of the latter from QCD is thus very important. The 
1/rriQ corrections suffer from uncertainties of different sorts, and a precise 
treatment is difficult. We have shown that the 1/mQ-corrections to elastic 
contributions can be estimated to a good extent, and, within the approxima- 
tions used here, for pseudoscalar mesons they turn out to be relatively small 
due to the fact that kinematic and hyperfine corrections conspire to cancel 
each other. For j3 = 3 GeV -1 the results are remarkably close to those ob- 
tained at leading order. To fully pin down the corrections to the elastic term 
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knowledge of the tuq dependence of Fy,3 and gy,s is required, and moreover, 
the difhcult-to-estimate tadpole terms should be included. Much more work 
is needed in order to estimate the importance of the disregarded inelastic 
terms. 

Our main aim was to conceptually clarify the analysis of electromagnetic 
splittings, and therefore, we refrained from making definite numerical claims, 
which are sensitive to the model chosen for the elastic form factors and other 
assumptions. Our estimates for the relative sizes of the 1 / itlq corrections to 
A 7el M# should however be taken seriously. 
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CAPTIONS 



Table: Results for electromagnetic mass splittings of pseu- 
doscalar mesons in the elastic approximation as defined in the 
text. Bold type numbers refer to I/tuq corrected values. The 
results correspond to the choice A ~ M# in eq. (19), and the 
error shown in the last column results from a change in A by a 
factor of 2. For comparison the results in the niq — > oo limit 
are also shown. The numbers designated by boxes are the results 
of the VMD model @. 



Figure: The diagram shows the flow of the large space like mo- 
mentum q in the process H — > H' + X mediated by the light 
quark component of the electromagnetic current. 
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